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SUMMARY 


A general method developed for the analysis of inviscid hyper- 
sonic shock layers is discussed for application to the case of the 
shuttle vehicle at high (65°) angle of attack. The associated extensive 
subsonic flow region caused convergence difficulties whose resolution 
is discussed. It is required that the solution be smoother than 
anticipated. 


INTRODUCTION 

The supply of solutions for the inviscid flow in a shock layer, such 
as that about a blunted body moving at hypersonic speeds, is very large 
and it might appear that another such method would be, to say the least, 
superfluous. However the goals of this work are special. The aim is to 
find an efficient and economical solution for the shock layer on the wind- 
ward side of the shuttle vehicle at 65° angle of attack. An exact analysis 
is not sought but only a reliable one. Thus the precise time -dependent 
methods or even the method of characteristics are rejected because of 
the required computing time. On the other hand, Newtonian methods are 
simple but do not provide the desired detail in the shock layer. 

Some time ago Maslen (ref. 1) proposed a method for axisymmetric 
flow based on the use of a Mises transformation coupled with a simple 
approximate integral of the lateral (to the shock) momentum equations. A 
number of solutions were given and the results were indeed accurate and 
very easy to obtain. 

The simplicity of the method led several others to apply the method to 
elaborations of the original case. Jackson (ref. 2) considered viscous 
effects and also used the method to solve the complete equations of motion 
iteratively. Perini and Melnik (ref. 3) solved nonequilibrium flow, while 
Olstad (ref. 4) applied the procedure to radiating flow and to the massive 
blowing problem associated with extreme heating rates. Recently, Grose 
(ref. 5) has presented detailed solutions for nonequilibrium hypersonic 
flows for a variety of cases of interest in planetary entry. 



Subsequently, Maslen (ref. 6) returned to the original inviscid 
equilibrium flow case, but considered general three-dimensional 
geometry. None of the elaborations mentioned above were examined. 
However the basic analysis was presented for a general smooth 
three-dimensional geometry. The present study is an extension of 
that analysis to the specific case cited above, the windward surface 
of the shuttle vehicle at high angle of attack. 
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SYMBOLS 

shock curvature functions, eq. (35) 
constants, eq. (36) 

auxiliary functions, eq. (5) or eq. (2Z) 
speed of sound behind shock 

enthalpy 

differential operator, eq. (7) 
free stream Mach number 

pressure 

pressure functions, eqs. (19), (20) 
entropy 

velocity components, cylindrical coordinates 

velocity components, defined in eqs. (2) 

cylindrical coordinates, fig. 1 

cylindrical coordinates (wind axes) 

shock layer coordinate, fig. 1, and eq. (1) 

angle of attack 

auxiliary function, eq. (28) 

ratio of specific heats 

auxiliary function, eq, (5) 

shock layer coordinate, fig. 1 and eq. (1) 
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e 


"Mises" type coordinate, eq, (6) 


X 

V 

P 

cr 

T 

Subscripts: 

s 

oo 


shock layer coordinate, fig, 1 and eq. (1) 

-I / j_ \ 

Tan ' r 3 0 fig. 1 
Tan ^ ( 9r/3z) , fig. 1 
flow angle, eq. (2) 
density 

"Mises" type coordinate, eq. (6) 

Tan ^ (Tan v Cos X. ), fig. 1 
stream function, eq. (3) 

stream function, eq. (3), also "Mises" type coordinate, 
eq. (6) 


shock value 
free stream value 


X. r, w 

z, ri, e 
O' . 4^ 


partial derivative 


GENERAL ANALYSIS 


The equation to be solved are the inviscid Euler equations for 
steady, compressible, isoenergetic flow. The development of the 
appropriate approximate system actually solved is discussed in reference 
6 and the following is a summary of the development in that report. 

Consider a set of coordinates (z, t), 0) related to the flow geometry 
(Fig. 1). We base them on the shock where z, 0 are the axial and azimuthal 
coordinates of a point on the shock while q is the (normal) distance from 
shock to field point. These coordinates are related to cylindrical ones 
(x, r, o)) by 
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( 1 ) 


X = 2 + T| Sin T 

r2 = + T]2 Cos^ t . 2Tlr Cos f Cos X 

s s 

r Sin Cl) -6) = Tl Cos t Sin X 

where equation of the shock and v > X, and 

T are angles defined in the symbol list._ Now introduce three mutually 
perpendicular velocity components, U, V Cos V Sin | where U is normal 
to the shock and i is otherwise arbitrarily chosen. Then 

U = Sin ^ “ ‘'^1 Cos T Cos (lO-S+X) + Wj^ Cos r Sin (to-9+X) 

V Cos ^ Cos T + Vj^ Sin t Cos (u>- 9+X) - Sin t Sin (cb-0+X) (2) 

V Sin ^ = “ Vj^ Sin (od- 0+X) - Wj^ Cos (u)-0+X) 

The Euler equations can be solved by the introduction of a pair of stream 
functions, ij; and <(> , such that 


pU = 7i|f x^cp 


(3) 


or 




pV{DCos^ + E Sin " ^'Tl '^2 

pV(ACos| + B Sin e) - r ' 
where A, B, D, E and A are geometric factors given by 


(4) 


A =“ [tg + Tl Cos t Cos X (Xg“l)]/r Cos X 
B = H T. 


D - Tan t Tan X + TlX Cos t 

' z 

E = (1 + T)t Cos t )/Cos t 
z 

A = AE-BU 


(5) 
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Now introduce a final coordinate system. This is a sort of Mises 
transformation in which we interchange *Jj and t| as dependent and 
independent variables. The specific use will clearly depend on the 
way in which i|/ and <r {eq. (3)) are subsequently chosen. Let the 
new independent variables be ij' > 

where 

^ = z 

0 - = e (6) 

ifj = z, 0) 


Further, define the operator L by 

L . -(D Cose + E Sine )±_ + (A Cos e + B Sin | ) 

8 £T 


8 

8^ 


(7) 


Then equations (4) yield, noting (5) 

_L (<p) = 0 

VL (ti) = U(AE-BD) 

pVn = —2^ = i 

^ ^4; r(A Cos ^ + B Sin^ ) D Cos! + E Sin § 

The energy equation becomes 

L(S) = 0 

while the isoenergetic condition is 

1 2 - 2 

h + (U + V ) = constant 


( 8 ) 

(9) 

( 10 ) 


( 11 ) 


( 12 ) 


The momentum equations become 


- ^ TT = pV [L(U) 


- V Cos ^ L(t) + V Sin ^ Cos t L(\)-G) ] 


(13) 
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( 14 ) 


- [A(P^ - DC/M) - “ Da/M)3 

= pV [ L(V Cos C) + UL( t) - V Sin C Sin t L(X-6)] 

- [B(P^ - DC/M) - f (15) 

= pV [ L(V Sin ^ ) + (V Cos g Sin t - U Cos t) L(\-0)] 


These last three equations can be combined to give the Bernoulli relation 


L(P) + I L(U^ + V ) = 0 


(16) 


The equations to be solved then consist of (8), (9), (10), (H), (12), any 
two of (13) -(16) and a state relation. The dependent variables are the 
three velocity components (U, V, i), the thermodynamic variables (P, 
p, S, h), one physical coordinate (r|), and one stream function (<P). 


For later use, 
some labor 


we will need the values of P 


4^ 


at the shock. 


After 


pr| 


± 


_ 2 

= V Cos T Cos X.[Cos C L.(t) 

1 


- Sin C Cos T Li(\- 0)] 

.2 


dU , U dP 

• 'T ■ ^ ^ 


U dSin 0 pU dSin p 
oo P ^ oo _J1 

1 - U^/c^ 


(17) 


2 r 

-2L , 

Poo 


Cos T Cos % \ Sin 
Cos T -CI-Xq) 


T Sin 1 Igl )l - p^/p) 

ITT?;7| 


where c is the speed of sound behind the shock. Note that all the 
derivatives on the right side lie along the shock since, clearly, £, and cr 
in the operator L as used here can be replaced by z, 0. 

As discussed in reference 6, we now approximate the momentum 
equation normal to the shock (eq. (13)) by its value at the shock, eq. (17), 
setting the terms in curly brackets equal to unity and assuming that the 
second line varies linearly with ijj across the shock. In that case one 
obtains 
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(18) 


P(KC,a) = Pg(C,a) + Pi<C,cj)('|f - 'Pg) - P2<C,a) 



where 


P. = 






V Cos T Cos \ 


[ Cos ^ L(T)-Sin ^ Cos x L,(\-G) ] j- 


and 

^2” [ Cos T T^-Cos T Cos \ (1-X.0)/ r- Sinx Sin \ Tq/ 


Finally, eq. (12) can be closely approxinnated as 


h 



constant 


while equations (7), (8) and (11) are applied, but setting xi = 0 in A, 
D, E. Thus, we have 

A = 1/Cos \ 

B := 0 


D Tan x Tan X. 

E = 1/ Cos X 

Then neglecting u and r\, the two momentum equations (14) and (15) 
be combined to yield 


pV^ [l<D + Sin T L(X - 9)1 


A Sin I P^ + (E Cos |-D Sin §) 


!=i 

r 


Finally, from eq. (7), the equation for the streamlines is 


^ ^ -(DV + EWl 
dQ r(AV + BW) 


(19) 

( 20 ) 

( 21 ) 

B. 

( 22 ) 

can 

(23) 

f24) 
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while along the streamlines^ eq, (23) becomes 


d Sin_i ^ 1 /p _ Tan T Sin X p \ Cos g Cos X 

dC I r V rCosr 


Cos X TanT (\g-i) 

+ Sin I - 




r Cos T 


( 25 ) 


Equations (18), (21), (8), (11) and (25) plus a state ^equation are a 
complete system of equations for P, p, h, S, c|), V and ^ . The only 
differential operator appearing is L (eq. (7)) which is the ordinary 
derivative along a streamline. 

After the system is solved, the solution is completed by returning 
to the physical (z, 0, r|) space via eq. (10). Experience in the axisymmetric 
case has shown that it is preferable to evaluate eq. (10) exactly without 
making the small Tj approximation of eqs. (22). 


Boundary Conditions 
In the free stream, we set 


cp = ou = Tan“^ r Sing) Co s q> - x Sin q> 
1 r Cos 0 ) 


(26) 


PQO^QO ^1 P 00^00 


2 2 
t (r Coso)) + (r Sin oo Cos a “ x Sin a) 1 


= Pqo^oq <» ] 

2 1 Cos cp J 


(27) 


where x^, r\ , cylindrical coordinates aligned on wind axis. Equations 

Equations (26) and (27) are consistent with eq. (3) in the uniform free 
stream and, for axial symmetry, is simply the usual Stokes stream function. 
The shock is taken to be Rankine-Hugoniot. The quantitites pU, V, 
and P p ^ continuous across it. Thus, immediately behind the shock 
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= Sin T Cos a - Cos t Sin O! Sin(0 -X) = Sin 


(28) 


pu/ Pqo ^QO 

V Cos £ /U = Cos T CosO: + Sin t Sin (X Sin (6-X) 
^ ' oo 

VSin£/U = -Sin OH Cos ( 6-X) 

r Sin 9 Cosof - z Sin a 
s 

r Cos 9 
8 



(29) 

(30) 

(31) 


’ ""00 ’ 


Cos a Sin p) - r (Cos t Cos X + Sin 0 Sin or Sin p)} 

(32) 


where Sin p is defined by eq. (28). Clearly, p is the complement of the 
angle between the wind direction and the normal to the shock so that the 
Mach number normal to the shock is M Sin p. The Rankine Hugoniot 
relations yield 


P/Pqq = 1 + Y^Sin^ P (1 - Poo^P^ (shock) 


(33) 


where, for a perfect gas 


2 ^ ^^00 ^ (Moo (shock) 


(34) 


At the body surface, we would like to have 4* = 0. This requires 
(see eqs. (22) and (26) that the origin of coordinates be chosen so that 
= 0 occur on the stagnation streamline. If this selection can be made, 
there is no further worry about body surface conditions. 
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METHOD OF SOLUTION 


The solatioa is found in an inverse manner. The shock 
curvature is assumed and is twice integrated to find the shock angle 
and radius. From that, the pressure can be found (eq. (18)). Then 
on a specified streamline, the entropy is constant at its value on 
crossing the shock. The state equation and isoenergetic condition 
(eq. (21)) yield the density, enthalpy and total velocity (V). Then 
eqs. (24) and (25) can be solved for the streamline position and flow 
angle. That completes the solution in the "Mises Space” ( ^ , cr , iji). 

We return to physical space, finding ti, the distance normal to the 
shock, from the first of eqs. (10). As i(j = 0 defines the body, the 
computed value of p at 4j= 0 must be compared to the geometric distance 
and iterated as necessary. A somewhat more complete description of the 
program is given in Appendix A. 

Implementing the procedure described requires two distinct parts. 
Once a satisfactory starting solution has been found, one can march 
downstream. Extrapolate the shock curvature to a new station and pro- 
ceed exactly as described above. 

For the starting solution, a different process is used. An initial 
distribution of the shock curvature is assumed over a region snrroianding 
the nose. Preferably this is the entire subsonic region (this point will 
be discussed further later on). This is taken in the form 

2 

curvature = a(cr ) + zb(<r ) + z c(<r) (35) 

where the most general form for a(cr) is 


a(cr ) = a 

1 + b Cos^o- 


(36) 


Then the solution is constructed for a set of z- stations (typically 8). 
The solutions at two intermediate stations (say 2 and 5) in two mutually 
perpendicular (in cr ) planes are used to find the values of a, b and the 
shock standoff distance at the nose. Then the values of b(cr) and c(o* ) 
are found iteratively by matching the body position and slope at the last 
(8th ) station. Once this starting solution converges, the downstream 
solution can proceed straightforwardly based on extrapolation of the shock 
curvature as discussed earlier. 


All of the elements described are straightforward and proceed without 
difficulty except for the iteration involved in determining the initial curva- 
ture distribution elements b(cr), c(cr) in eq. (35). 
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Originally, the initial distribution was made with c(cr) 3 e 0 in 
eq. (35). In that case only the body position was matched (at station 8), 

For some problems, this appeared to be acceptable and the solution 
could continue downstream without especial difficulty. It was observed 
however that the resulting surface values of, say, pressure or Mach 
number were not smooth but had a bump just downstream of the starting 
solution. The same behavior had been observed in a code developed 
for axisymmetric flow. This corresponds to the fact that the body, 
slope is poorly matched at that point and the resulting corner’^ generates 
some mild havoc. 

Unfortunately, there was a very large snake in the grass. The 
problem posed in the present study is not typical of the cases for which 
the bumpy solution could proceed. On the windward surface of the 
shuttle vehicle at 65*^ angle of attack, there is a very extensive sonic 
or near sonic region. It is not really practical to have the initial solution 
cover this whole region. Nevertheless, the initial solution can still be 
found. However one must then march downstream toward the sonic 
point'^. The solution so found is divergent at a critical point near the 
limiting characteristic separating the sub and supersonic regions. The 
saddle point behavior is familiar from the method of integral relations 
to which the present method bears considerable resemblance. 

Unless the initial solution is a good one, the solution resists the 
approach to the critical point with great enthusiasm. Accordingly it 
became clear that a smooth solution would be needed. For that reason, 
the third term (c(cr ) in eq. (35)) was introduced. Then both the body 
position and slope can be matched at the end of the initial region. For 
axisymmetric flow, this added term is no great complication and its 
introduction serves to give good results even when a subsonic region lies 
downstream. In the present case, the selection of a reasonable iteration 
procedure becomes of paramount importance. Clearly we cannot perform 
indefinitely many iterations and hope to have an efficient solution of the 
problem. 

A number of procedures have been tried. First several brute force 
calculations were made to clarify the expected behavior of b(cr) and c(cr) 
when good convergence is obtained on the body position and slope. Next 
some limited runs were made to get an idea of the precision required in 
satisfying these matching conditions if the subsequent solution were to 
proceed without difficulty. These calculations were done for an hemisphere 
cylinder at Mqo ~ 10*0 and y- 1.4 for a perfect gas. Angles of attack of 
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0, 40 and 50 degrees have been run. The 0 degree case can be 
compared with similar results from an axisymmetric code. The 
results agree to within about l/ 2% for shock shape and surface 
pressure for the two programs and are in excellent agreement 
with other calculations. The case serves largely as a check on 
the overall operation of the present general program. At 50° angle 
of attack, the windward flow is near sonic for some distance and is, 
in that sense, typical of the desired shuttle condition. While great 
improvement was found in calculations using the complete eq. (35), 
convergence of the iterative process was not good enough. Time 
did not permit a satisfactory exploration of this vital question. To 
reduce the sensitivity to sonic behavior somewhat, an angle of attack 
of 40° was also run. While the iteration scheme which was used con- 
tained a strong element of trial and error, the results are gratifying 
in that reasonable convergence was obtained and it became clear that 
the solution could continue downstream. 

An extended discussion of some of the problems associated with operation 
is given in Appendix B. Particular attention is directed to the effects of 
a sonic region. 
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CONCLUDING REMARKS 


The present study has shown that the analysis of three dimensional 
blunt body flows can be carried out by the method based on approximate 
integration of the momentum equation across the shock layer. However 
it is clear that considerable care and more study is required to provide 
a suitably rapid convergence of the required iterations if an efficient procedure 
is to be attained. The exploration of such iterations thus far carried 
out indicates that success is to be expected. 
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APPENDIX A 


The computer program consists o£ a main program plus eleven 
subroutines. These are: 


NONAX 

Main logic program controlling the sequence 
of use of the subroutines. 

CHEM 

Does all chemistry calculations for a perfect 
gas. Could be replaced by real gas 

START 

Provides initial values at the nose as well as 
starting estimates of the flow near it. 

TEST 

Routine to test (and provide an iteration to 
improve) the starting solution. 

SHOKQ(I) 

Computes quantities at the shock for specified 
radius, slope and curvature distribution. 

INTEGG(I) 

A third order Runge-Kutta routine to integrate 
the equation for the streamline (eqs 24 and 25) 
(one step for each streamline). 

AITKEN 

Interpolation routine for use in INTEGG 

TRAN(I) 

Transfer the results of INTEGG from the o- , i|< 
coordinates to z, 0, t) ones and integrates 
the first of eqs (10) to find ti. 

CROS(I) 

Tests the position of the body from TRAN 
against the actual geometric coordinates and 
provides an iteration scheme. 

EXTRAP 

Extrapolates the solution one step downstream 
when the results in CROS converge. 

EXTEND 

Routine to extrapolate a small distance into the 
separated leeside region. 

SMOOTH 

A least squares smoothing routine to smooth the 
shock estimates so that two lateral derivatives 
can be found. 
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A listing o£ these routines is not presented (although available) 
because their present form contains a confusing number of logic elements 
associated with various iteration procedures which will have no part in 
a final program. However the general sequence is to make an initial 
estimate of the shock slope and position (START), use that to find shock 
values for 5 stations (SHOKQ), integrate the streanaline equations (INTEGG) 
and transform to physical space (TRAN) and test the result (TEST, CROS) 
to determine an acceptable nose curvature distribution (eq (36)). Then 
one goes through INTEGG, and TRAN to station 8, uses CROS (8) to 
test the body position and slope at station 8. Once that converges the 
solution is extrapolated (EXTRAP) and one goes thru SHOKQ, INTEGG, 
TRAN and CROS, interating as necessary. 
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Appendix B . 


This section discusses some of the details of the agonies 
connected with getting the program to run and is in the nature of a 
status report on its present imperfect condition. 

First of all, all of the subroutines which compute the consequences 
of an assumed shock curvature distribution (starting routine) or 
extrapolation (downstream continuation) are satisfactory. These 
comprise the construction of the shock slope and radius, evaluation of 
quantities behind the shock (SHOKQ) and hence the pressure distribution, 
integration (INTEGG) of the equations for the streamline and flow direction 
(Eqs. 24 and 25), transformation of these results to ^ = const lines and 
integration of eq. (10) to find )r^ and thus return to physical space (TRAN) 
and then to compare (CROS) the body coordinates thus generated with the 
actual ones found geometrically. 

The operating status and attendant problems can be seen with the 
aid of figs. 2 and 3. These give the shock shape, and surface pressure and 
Mach number distribution for an hemisphere cylinder solved for a perfect 
gas at M = 0 and ^ - 1. 4. In figure 2 the angle of attack is 40°, In this 
case the starting solution applies to station 8 and yields a solution to that 
point which is entirely supersonic and which involves only the spherical nose 
of the body so that the starting solution remains axisymmetric with respect 
to the wind directions. In that case the starting solution is one in which 
the body position is matched at station 8 correct to • 004 (out of .18) 
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while the slope is 3% below the geometric value. The subsequent 
downstream calculation proceeds without any difficulty about convergence 
(a maximum of 3 iterations per station). The calculation was halted 
arbitrarily, not because of any operational problem. For future discussion, 
the results on the windward streamline should be noted. The shock surface 
pressure and surface Mach number vary smoothly and, importantly, 
the surface Mach number is about 1. 2. Off the surface (not shown) the 
values are larger; increasing monotonically toward the shock. 

Now consider figure 3, which shows corresponding results for 

= 50^, In this case, the starting solution covers a small portion of 
the cylindrical afterbody. As before, at Station 8, the solution has position 
correct to . 005. (out of . 18). 

The error in slope is about 7% on the windward streamline. There 
was no problem with the solution to this point. On marching downstream 
the solution (fig. 3) has no unusual features over any but the windward part 
of the body. Except for that region, the results are as one would expect as 
one proceeds from O^: = 40^^ (fig. 2) to 0^=50°. 

However, on the wind side there is trouble. The solution found is as 
indicated, with, in particular, the Mach number on the windward stream- 
line becoming subsonic and falling to relatively low values. Under these 
circumstances the error criterion on (body position) was never satisfied 

at stations 11 or 12 (for Q = 90, 85). The program accepted this solution 
when no improvement occurred after several iterations. 

The principal difference between the = 40 and 50 cases is the 
indication in the latter case of a local near sonic region. It has been shown 
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in other analysis (ref. 7 and this author ’ s unpublished work) that the 

present method is much like that of integral relations in that the 

problem is singular near a ’’sonic” region. This arises because (eq. 10), 

the changes in the value of with changes in the shock geometry have 

an effective influence coefficient roughly proportional to a weighted 
2 

integral of (M -1) across the shock layer. In two dimensional or axi- 
symmetric flow this behavior is very straightforward and clear; in three 
dimensions the consequences are rather obscure. In any case, there is 
danger. 

This describes the status and raises some specific questions which 
can be examined in turn. These include; 

1. Is there a fundamental problem due to violation of required 
characteristic propagation directions? 

2. Is the solution sensitive to the position of Station 8? 

3. Must the slope in the starting solution (Station 8) be matched 
much more accurately than at present? 

4. Is there a local ’’sonic” region which can be "jumped”? 

5. Is the solution, in fact, correct? 

6. What is the effective way to iterate b, c (eq. 35)? 

7. Do we require a body with continuous curvature? 

At the present time, none of these questions can be resolved entirely, 
but some indications are available. The questions about propagation of shock 
effects toward the body is not a local error. In the present case, one really 
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has the shock affect propagated normal to the shock and hence nearly 


normal to the body which, from the viewpoint of characteristics is 
most valid if M is near 1. Actually the local shock effect propagated is 
rather diffuse since we specify the shock by its curvature; that is, we in 
effect require it to be twice differentiable. There seems to be no unusual 
error due to this behavior, and if anything the near sonic regions are most 
realistically portrayed in this respect. 

The sensitivity of the results to the specified position of Station 8 
(Question 2) has not really been examined. A few runs have been made (at 
low angle of attack) with that station moved aft about 10%. The results were 
unchanged. No other tests have been run to date. It should be noted however 
that in the c?(= 50 case (fig. 3), the entire flow at Station 8 was supersonic. 

It was not just the "average” value. 

The question (#3) of the effect of slope errors is less clear. In the 
particular run given in figure 3 the error in position and slope is as sketched 
below. In particular, the two errors are in a direction to tend to compensate 
each other on going to Station 9. This would then require recompression in 
going to Station 10. 



£m>rs ccfSkMws ? 
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This brings as to question 4. If there is a real ^*near sonid^ 
region, can one jump over it or is it too extensive for a meaningful 
jump? This question is, in the end, related to question (5); what is 
the real answer? No other results appear to be available* However 
it seems pretty unlikely that the result of figure 3 is correct at least 
so far as it shows a large drop of the Mach number (to * 80 at Station 12) 
on the windward streamline. 

Question 6 was discussed earlier in the main text. It may be 
that one can iterate (b) generally but only require that (c) match formally 
at a few points. This remains open. 

Finally, the body tested has discontinuous curvature at the juncture 


of the spherical nose and ^cylindrical afterbody. 


Such a jump is a frequent 


source of difficulty with flow fields and is usually resolved by a fillet at 


the juncture forcing continuous curvature. In the present analysis this 


may not be a problem because of the relatively large step size (about 10% 


of the nose radius) in that region. It should be noted however, that the 


surface pressure irregularities in the computed examples (figs. (2) and 
(3)) occur at about the juncture of the nose and afterbody. 
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(c) Shock Shape 


Shock Shape and Surface Pressure and Mach Number for Hemisphere 
Cylinder (0« = 40°). 









